Quarkonium-Hadron Interactions in QCD 



Yongseok Oh,0 Sungsik Kimj^ and Su Houng Lee 



(N 
O 

o 

Ph. 
< 



(N 
> 
(N 
m 



^ ■ 

Oh, 
I 

Oh! 



X 



Institute of Physics and Applied Physics, Department of Physics, Yonsei University, Seoul 120-749, Korea 

The leading order non-relativistic QCD result for the quarkonium-hadron interaction cross section 
is derived by using the QCD factorization formula combined with the Bethe-Salpeter amplitude for 
the heavy quark bound state. Our result shows the equivalence of this approach with the operator 
product expansion even for the bound state scattering to leading order in QCD. We also apply 
the formula to the charmonium system and investigate the relativistic correction coming from the 
phase space of the reaction. We find that although JZ-i/i-gluon cross-section gets a large relativistic 
correction at higher energies, the correction becomes small in the J/i/)-hadron cross section due to 
the increase of the gluon distribution at smaller x inside the hadron. 

PACS numbers: 25.75.-q, 12.38.-t, 13.75.-n, 13.85.Lg 



More than two decades ago, Peskin and Bhanot 
and Peskin showed that one could apply perturbative 
Quantum Chromodynamics (pQCD) to calculate the in- 
teractions between the bound state of heavy quarks and 
the light hadrons. Such calculation was feasible because 
in the large quark mass limit one could consistently ob- 
tain the leading order operator product expansion (OPE) 
of the correlation function between two heavy meson cur- 
rents in the light hadron state. The justification of such 
calculation stems from the fact that there exist two rel- 
evant scales in the bound state Q in the large quark 
mass limit, namely the binding energy, which becomes 
Coulomb-like and scales as mg^, and the typical mo- 
mentum scale of the bound state, which scales like mg^, 
where m is the heavy quark mass and g the quark-gluon 
coupling constant. Hence, taking the separation scale of 
the OPE to be the binding energy, it is consistent to take 
into account the bound state property, which is generated 
by the typical momentum scale of the bound state, into 
the process-dependent Wilson coefficient. 

The calculation is challenging in itself, and has at- 
tracted recent interests because when applied to the J/ip 
system |^, the strength of the J/i/i-hadron inelastic 
scattering cross section determines the amount of the J/^p 
suppression |^ in heavy ion collisions due to hadronic fi- 
nal state interactions Q . Moreover, other existing model 
calculations for (Jj/^-}n based on meson exchange mod- 
els [H, quark exchange models [||, or QCD sum rules 
Ip^ give very different energy dependence and magni- 
tude near threshold, which shows the importance of the 
nonperturbative effects. Therefore, more careful analyses 
of each calculation are necessary, where the limitations 
and expected corrections would be explicitly spelled out. 

In this paper, we will derive the leading order pQCD 
result using the QCD factorization theorem. The result 
obtained by Peskin and Bhanot was derived within the 
OPE H, |2| . The equivalence between the OPE and the 
factorization approach is well established in the deep in- 
elastic scattering and Drell-Yan processes . As will be 
shown here, the equivalence between the two approaches 



are also true for bound state scattering to leading order in 
QCD. Similar approach has been used by one of us (SHL) 
to estimate the dissociation cross section of the J /ip at 
finite temperature The factorization formula also 

provides a manageable starting point to calculate higher 
twist gluonic effects which should be nontrivial 

for the J/V'-hadron scattering. 

We refer to a bound state of heavy quark and its own 
antiquark as $. According to the factorization formula, 
the total scattering cross section of $ with a hadron h 
can be written as 0, p3 
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with V = p ■ q/M>j,, where p (q) is the momentum of the 
hadron ($) and A/$ is the $ mass. Here, cr^g is the 
perturbative <i>-gluon scattering cross section and g{x) is 
the leading twist gluon distribution function within the 
hadron. The separation scale is taken to be the binding 
energy of the bound state. Hence, the bound state prop- 
erties have to be taken into account in cr$g. This can 
be accomplished by introducing the Bethe-Salpeter (BS) 
amplitude r(pi, — P2), which satisfies [l^ 



r/.(pi 



-P2) = iCcolor 

X r^(pi + k, 



g^V{k)YMPi + k) 



P2 -f A:)A(-p2 + fc)7i^, (2) 



where Ccoior — {N^ — l)/(2-/Vc) with the number of color 
Nc- The trivial color indices have been suppressed. Fig- 
ure |l] shows the diagrammatic representation of the BS 





FIG. 1: The Bethe-Salpeter equation for the quarkonium 
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equation, where pi (— P2) is the four-momentum of the 
heavy quark (anti-quark) with A(p) = l/^^—m + ie) and 
V{k) = -l/ik^ + ie). 

We introduce 6,,, which is defined as 



Nc_ 



1/2 



2tt 



(3) 

where q = pi +P2 and p ~ {pi —p2)/2. Wc work in the $ 
rest frame and introduce the binding energy e, such that 
qo = M$ = 2m -I- e (e < 0) . Then, in the non-relativistic 
hmit, Eq. (p|) reduces to 



e — p^/to V M, 



Nc 1 + 70 



xr,[-+p,--+p^ 2 



1 - 70 



(4) 



where po = e/2 — p^/(2m). In this hmit, the BS amph- 
tude reduces to 



r^{q/2+p,-q/2 + p) 



-70 . 1 



70 



(5) 



and the corresponding BS equation becomes the non- 
relativistic Schrodinger equation for the Coulombic 
bound state, 
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y(k)^(p + k), 



(6) 

so that V''(p) is the normahzed wave function for the 
bound state. For the pseudo scalar bound state, we re- 
place 7i5^i to 75 in Eq. (||), but obtain the same spatial 
wave function. (Some discussions on the realistic poten- 
tial can be found, for example, in Ref. [0.) It is straight- 
forward to verify that our formalism gives the well-known 
result for the $ e+e~ decay width, r(<I> — > e+e^) — 
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\^p{r = 0)1 , where ac 



7(47r) and Qh is 



the heavy quark charge in the unit of e. 

With the BS amplitude defined as Eq. (||) , we now ob- 
tain the $-gluon scattering amplitude with the processes 



depicted in Fig. ||. The scattering matrix elements can 
be written a,s M = £p($)A^^''e,y(.g), where e^($) and 
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FIG. 2: Scattering processes of $ and gluon. 

£1/(5) are the polarization vectors of the $ and the gluon, 
respectively. We first compute M.i'^, which comes from 
Figs. |(a,b), 

= -.g^l(pi)[r^(pi,-<z+Pi)A(-g + pi)7'^ 
+ 7''A(g - P2)T^{q - P2, -P2)]T"i;(p2), (7) 

where u{p) and w(p) are the Dirac spinors of the quark 
and antiquark, respectively, a is the color index of 
the incoming gluon, and is the color matrix with 
Tr(r"T'') = 5"V2- To obtain the leading order non- 
relativistic result, we have to discuss the counting scheme 
in the rest frame of First the binding energy = 
m\Nc9'^ is O(TOg^). Combined with the energy 
conservation q + k — pi + p2 in the non-relativistic limit 
this implies the following counting, 

|pi|~|p2|^0(mg2), ^\k\^0{mg'). (8) 

Then the leading order result of Eq. (^ reads 
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(9) 



Figure l|(c) gives 



Mr = g'fabcT'T''u{p,) J ^V{l-p,)V{l-p, + k)^fjA{l)r'^{l,l-q)A{l-q)-f^vip2) 

X [g"''{2l - 2pi + kY + g^-'ipi ~ I + k)" - g'^'il - pi + 2kf] . (10) 

I 



The color factor of this diagram is fabcT^T^ = -fiVcT". 
There is one more diagram, where the external gluon leg 
is attached to the internal quark line within the inter- 



nal loop. However such a diagram carries the color fac- 
tor T^T'^T'' = -T°-/{2Nc) (color adjoint potential) and 
hence is suppressed by compared to Fig. ||(c) |jl|. 
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Using the BS equation and the counting (H), we get leading order result for the scattering amplitude as 



x^bi)^y^^r"t'b2). (11) xu(pi)l±^yi^r««fe). (12) 



Collecting Eqs. (^ and (|ll|) gives the (gauge invariant) The scattering cross section cr$g can be obtained from 
I 
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M$(fcoP? - m2) + m2(p0 - fcn) +P?(M| - 2m2)/2 Gin'^AUko 
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IXPdriiab (13) 



in the $ rest frame, where jA^p — 6{n^—i) 

EIA^P and 

the average is over the initial polarizations and the color 
of the gluon. Note that going from the second line to the 
third hne in Eq. (|l|), we use the non-relativistic approx- 
imation (|^). With the amplitude ( p^ ) we obtain 



|X|2 = 



vv(p)r 



(14) 



Making use of the ground state wave function of the 
Coulombic bound state, 



VV'is(k) 



5/2 



32^/¥■ 



(15) 



where e — — Eq and ao — IQi: j [g^ N cm) , we finally obtain 



cr$g(A) 



128g2 ^(A/eo- 1)3/2 



(16) 



with A = g • k/Mq,. Hence, we have derived anew the 
result for a$g of Ref . . The scattering cross section cr^/i 
can then be computed using the factorization formula (|^). 

To apply the same procedure to the scattering of $' 
(the 25* state), we use the wave function 



VV2s(k) 
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flok, (17) 



with the binding energy e = — eo/4. This then leads to 

(18) 



165^ ,iX/so-\r^'iX/eo~ir 



{X/eoY 



which is consistent with the result of Ref. ||^ . It should 
be noted, however, that the application of our formalism 
should be less reliable for the excited states becasue the 
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FIG. 3: Scattering cross sections o",j/^_p and cr^'{2S)-p- 
solid lines are obtained with g(x) of Ref. [^] (BP), while the 
dashed and dotted lines are with that of Ref. jis) (MRST) 
13 (GRV-LO), respectively. The cross sections with rel- 



The 



and 



ativistic corrections are given by dot-dashed lines with g(x) 
of Ref. Ili (GRV-LO*). 



wave function would be more affected by the confining 
part of the potential. 

We now plot the numerical results for the (T$p and cr$/p 
in Fig. ^. Similar plot has been already given in Refs. 
iQ, |j, but here we use the original parameters as given 
in Ref. [§: eo = 780 MeV and to = 1.95 GeV as deter- 
mined in Ref. Q by fitting the J/?/' and ip'{2S) masses 
to a Coulombic spectrum. Also given in Fig. || are the 
results with the relativistic correction (dot-dashed lines) 
that will be discussed below. For the gluon distribution 
function of the proton, we use the parameterization of 
Ref. @, g{x) = 0.5{r] + l){l-x)i/x with 77 = 5 (BP) as 
well as the MRST |l|] and the GRV-LO gluon distribu- 
tions |l^. The results show the dependence of the cross 
sections on the gluon distribution function. With the 
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distribution functions of Refs. ^ the cross sections 
increase with the center-of-mass energy ^/s. But g{x) of 
Ref. [|8| leads to the decrease of the cross sections at 
higher energies, which comes from the difference in g{x) 



at small x. Looking at the ratio cr ,, 



i/>'(2s)-p/^J/i/'-p 

this approach, one finds that it has a peak at ^/s « 4.2 
GeV. The peak value is about 2000 ~ 5000 but it quickly 
converges to 4 20 at higher energies depending on the 
form of g{x). In a quark exchange model the ratio be- 
tween the maxima of the cross sections was estimated to 
be around 6 Istl. 
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FIG. 4: 



The 



Scattering cross sections aj/^_g and <y^i(2S)-g- 
solid lines are obtained in the non-relativistic limit and the 
dashed lines are with the relativistic correction. 



Let us now consider the relativistic corrections to our 
results. The non-relativistic counting scheme (^) of this 
formalism implies that the formula ( p^ ) should not be 
used at high gluon energy. However, calculating rela- 
tivistic corrections in this formalism is intricately related 
to the higher corrections, which is a formidable future 
task. Here, we will only look at the simple but impor- 
tant relativistic correction coming from the relativistic 
calculation of the phase space integral in Eq. (|l3|). The 
difference after a full calculation is shown in Fig. ^ for 
cr$g , which shows larger suppression of the cross sections 
at higher energy. In order to see the relativistic correc- 
tions in cr$p we compute (T$p with the (J$g after the rela- 
tivistic correction using the GRV-LO gluon distribution 
function (the dot-dashed lines in Fig. ^). We find that 
this relativistic correction does not have much effect on 
the cross section cr$^. This follows because cT#g is highly 
peaked at small energy region and the gluon distribution 
function increases towards small x. Through Eq. (^, this 
implies that a^h, at all energies, are dominated by the 
low energy behavior of (T$g, which is reliably calculated 
in the non-relativistic limit. In fact, it has been observed 
in a potential model that the relativistic treatment of the 
bound state wave functions does not improve much the 
nonrelativistic description All these imply that the 
full relativistic corrections to the scattering amplitude 
(T$g would have little effects on cr^/i. 

The remaining important and interesting correction is 



the higher twist-effect, which contributes as A/e^ in our 
formalism. Since A is related to some hadronic scale of 
the hadron it should be of the order itself, hence 
should be non-negligible. The research in this direction 
awaits more experimental and theoretical works [^3[ Q . 
It is also important to extend this calculation to inves- 
tigate the absortion cross section for the x states. This 
will influence the amount of J /ij) production coming from 
the decay of the x's in a p-A or A- A reaction 



We are grateful to A. Hayashigaki and C.-Y. Wong for 
fruitful discussions. This work was supported in part by 
the Brain Korea 21 project of Korean Ministry of Edu- 
cation, KOSEF under Grant No. 1999-2-111-005-5, the 
Yonsei University Research Grant, and the Korean Min- 
istry of Education under Grant No. 2000-2-0689. 



skim@phya.yonsei.ac.kr ; Present 



Electronic address: yoh@phva.vonsei.ac.kr _ 
Electronic address: 
address: Department of Physics, Sungkyunkwan Univer- 
sity, Suwon 440-746, Korea 
Electronic address: 



[1] 
[2] 

[3] 

[4] 

[5] 
[6] 
[7] 
[8] 



suhoung@phya.yonsei.ac.kr 



[9] 



[10] 
[11] 

[12] 

[13] 

[14] 
[15] 

[16] 

[17] 
[18] 
[19] 
[20] 
[21] 



M. E. Peskin, Nucl. Phys. B156, 365 (1979). 
G. Bhanot and M. E. Peskin, Nucl. Phys. B156, 391 
(1979). 

T. Appelquist, M. Dine, and I. Muzinich, Phys. Rev. D 
17, 2074 (1978). 

D. Kharzeev and H. Satz, Phys. Lett. B 334, 155 (1994); 
Kharzeev et al, ibid. 389, 595 (1996). 
Arleo et al., Phys. Rev. D 65, 014005 (2002). 
Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 
Vogt, Phys. Rep. 310, 197 (1999). 
G. Matinyan and B. Miiller, Phys. Rev. C 58, 2994 

(1998); K. L. Haglin, ibid. 61, 031902 (2000); Z. Lin 
and C. M. Ko, ibid. 62, 034903 (2000); Y. Oh, T. Song, 
and Su H. Lee, ibid. 63, 034901 (2001); A. Sibirtsev, K. 
Tsushima, and A. W. Thomas, ibi d. 63. 044906 (2 001); 
W. Liu, C. M. Ko, and Z. W. Lin, ^lucLth/0107"058 . 

C. -Y. Wong, E. S. Swanson, and T. E. Barnes, Phys. Rev. 
C 62, 045201 (2000); 65, 014903 (2002); C.-Y. Wong, 
ibid. 65, 034902 (2002); private communications. 

F. S. Navarra et al, Phys. Lett. B 529, 87 (2002) 

B. Humpert and W. L. van Neerven, Phys. Lett. B 102, 

426 (1981); Phys. Rev. D 25, 2593 (1982). 

T. H. Hansson, Su H. Lee and I. Zahed, Phys. Rev. D37, 

2672 (1988). 

J. Bartels et al, hep-ph/9908411 (unpublished); J. Bar- 
tels and C. Bontus, Phys. Rev. D 61, 034009 (2000). 
S. Kim and Su H. Lee, Nucl. Phys. A679, 517 (2001). 
R. K. Ellis, W. Furmanski and R. Petronzio, Nucl. Phys. 
B207, 1 (1982); B212, 29 (1983). 

E. E. Salpeter and H. A. Bethe, Phys. Rev. 84, 
(1951). 

D. Ebert et al, Phys. Rev. D 62, 034014 (2000). 
A. D. Martin et al, Eur. Phys. J. C 4, 463 (1998). 
M. Gliick et al, Z. Phys. C 67, 433 (1995). 
W. Lucha et al, Phys. Rev. D 46, 1088 (1992). 
L. Gerland et al, Phys. Rev. Lett. 81, 762 (1998) 



1232 



